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The tunneling ionization of exotic atoms such as muonic hydrogen, muonium and positronium in a strong 
laser field of circular polarization is investigated taking into account the impact of the motion of the center of 
mass on the the tunneling ionization dynamics. The momentum partition between the ionization products is 
deduced. The elTect of the center of mass motion for the momentum distribution of the ionization components is 
determined. The effect scales with the ratio of the electron (muon) to the atomic core masses and is nonnegligible 
for exotic atoms, while being insignificant for common atoms. It is shown that the electron (muon) momentum 
shift during the under-the-barrier motion due to the magnetically induced Lorentz force has a significant impact 
on the momentum distribution of the atomic core and depends on the ratio of the electron to the atomic core 
masses. 

PACS numbers: 32.80.Rm,33.60.+q,36.10.Ee,36.10.Dr 


I. INTRODUCTION 

In a strong laser field ionization of an atom takes place by 
absorption of multiple laser photons HHH. The photons carry 
momentum which is distributed among the ionized electron and 
the atomic core (ion) after the interaction. The photon momen¬ 
tum transfer in ionization is a nondipole effect, theoretically 
described by a first relativistic correction to the nonrelativis- 
tic Hamiltonian. In classical terms the momentum transfer to 
the electron along the laser propagation direction arises due 
to the magnetically induced Lorentz force. At rather strong 
laser fields the ionization is in the tunneling regime, when 
the Keldysh parameter y is small, and the electron is released 
from the atom by means of tunneling through the potential 
barrier formed by the laser suppressed atomic potential. In 
the tunneling regime the Lorentz force induces the momentum 
transfer to the electron along the laser propagation direction 
either when the electron moves in continuum after releasing 
from the atom, or during the tunneling step of the ionization. 
The first effect is characterized by the the relativistic field pa¬ 
rameter f s eEdlmeCU), where the Eq and u are the laser field 
amplitude and angular frequency, respectively, —e and nie are 
the electron charge and mass, respectively, and c is the speed 
of light. The momentum transfer along the laser propagation 
direction in continuum is ~ t/p/c, where Up - 12 

is the electron ponderomotive energy in the laser field. It is 
responsible for shifting of the angular distribution of photoelec¬ 
trons from the laser polarization direction into the propagation 
direction in relativistically strong laser fields when ^ ~ lEIll, 
for the suppression of the electron rescattering with the ion 
151 and, consequently, suppression of nonsequential double 
ionization mm and high-order harmonic generation MB- 

The relativistic theory for the under-the-barrier dynamics 
demonstrated Da that the magnetically induced Lorentz force 
brings about a momentum shift along the laser propagation 
direction also during the under-the-barrier motion. It has a 
characteristic value of /p/c, which can be estimated from Apj. ~ 
eiValc)BoTii ~ Ip/c, with the atomic velocity Va = yjllpime. 


the ionization potential Ip, the laser magnetic field Bq = Eq, 
the Keldysh time tk — jlu, and the Keldysh parameter y - 
y[hJ^P- 

At nonrelativistic intensities the photon momentum transfer 
is rather small, nevertheless, it has been measured in a recent 
remarkable experiment IfTTI . using nonrelativistic laser intensi¬ 
ties of 2 - 10 X 10*"^ W/cm^ and detecting electron momenta 
of an order of 10“^ a.u.. The experimental results DD in a 
circularly polarized laser field indicated that after tunneling ion¬ 
ization the ion carries a momentum /^/c along the momentum 
of incoming photons, while the electron carries a momentum 
of t/p/c. This is in accordance with calculations of E3 and 
with the so-called simpleman model E3. The latter assumes 
that the electron appears at the ionization tunnel exit with a van¬ 
ishing momentum and then is driven solely by the laser field. 
The ion momentum in this case follows from the momentum 
conservation law. 

However, due to the Lorentz force effect during the under- 
the-barrier dynamics, the electron appears at the tunnel exit 
with a nonvanishing momentum along the laser propagation 
direction Da. The latter contributes to the asymptotic elec¬ 
tron momentum because of which the final ion momentum 
decreases with respect to the prediction of the simpleman 
model. Thus, the peak of the electron momenrnm distribu¬ 
tion along the laser propagation direction is at pe^ - Ipl3c (if 
the Coulomb field of the atomic core is neglected during the 
under-the-barrier motion), in contrast to the simpleman model 
vanishing prediction, and the ion should carry a momentum 
Pi^ - 21 p 13c Da El- Moreover, the ion momentum increases 
when the Coulomb field of the atomic core is accounted for 
in the near the over-the-barrier ionization regime, e.g., it is 
Pij = 0.8/p/c at EqIEu - 0.05, see Fig. 6 in ifTSl . where 
Ea - (2/p)^^^ is the atomic field strength. As it is noted in DI, 
these results are not inconsistent with the experiment of Ref. 
DD because of the large experimental error. The experiment 
on the momentum partition at ionization is most clearly shown 
in a circularly polarized laser field because in this case the rec¬ 
ollisions are avoided. In a linearly polarized field the Coulomb 
focusing because of recollisions modifies significantly the final 
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momentum distribution IfT^ and complicates the analyses HU. 
Recently, a numerical relativistic calculation of the electron 
momentum distribution in a linearly polarized laser field IfTSl 
indicated that the momentum shift in the laser propagation 
direction depends linearly on the laser intensity, which is in ac¬ 
cordance with the experiment CD. An additional momentum 
shift can arise also because of the, so-called, tunneling time 
GMD but this momentum shift is in the direction transverse 
to the laser propagation direction ifTSl I22l425]l and will not be 
discussed in the present paper. 

If the atomic degree of freedom is neglected during tunneling 
ionization of the electron, the momentum transfer to the ion 
can be simply deduced via the momentum conservation law, 
taking into account the electron momentum after the ionization. 
In treating the electron dynamics in strong field ionization, 
usually, the ion is assumed to be not moving. While it is a 
good approximation for the ionization of hydrogen and other 
common atoms because of smallness of the mass ratio of the 
electron to the ion, it is not the case for exotic atoms such 
as muonium [the bound state of an electron and antimuon, 
the mass ratio of an electron to muon is ^ 1/207], 

muonic hydrogen [the bound state of a muon and a proton, 
the mass ratio of a muon to proton is « 0.1126] and 

positronium [the bound state of the electron and positron, the 
mass ratio of an electron to a positron is me-lme+ = 1]. The 
hydrogen ionization problem taking into account the ion degree 
of freedom first was considered in |[T4l for derivation of the 
ion momentum, however, neglecting its impact on the electron 
momentum distribution. 

Note that the influence of the nuclear degrees of freedom on 
the electronic dynamics is well known in strong field molecular 
processes as, so-called, non-Born-Oppenheimer dynamics, see 
e.g. II26M1. 

In this paper we investigate the influence of the atomic core 
(ion) degree of freedom on the electron (muon) tunneling dy¬ 
namics in a strong laser field of a circular polarization in the 
case of hydrogen, muonium, muonic hydrogen, and positro¬ 
nium atoms. The impact of this effect on the photo-electron 
(-muon) momentum distribution and, respectively, on the par¬ 
tition of the photon momentum transfer between the electron 
(muon) and the atomic core (proton, antimuon, positron) are 
studied. 

We label by “electron” the negative component of the ex¬ 
otic atom, which is the electron in the case of muonium and 
positronium, and muon in the case of muonic hydrogen. The 
label “ion” is employed for the positive component of the ex¬ 
otic atom, which is the proton in the case of muonic hydrogen, 
antimuon in the case of muonium, and the positron in the case 
of positronium. 


II. THE SIMPLEMAN MODEL AND THE 
ENERGY-MOMENTUM CONSERVATION 

In this section we derive the momentum partition between 
the electron and ion during tunneling ionization, assuming 
that the photoelectron dynamics follows the simpleman model 
ina. The information provided by energy and momenrnm 
conservation laws is used to deduce the ion momentum p, and 
the number of absorbed laser photons n corresponding to a 
certain photoelectron momentum in the tunneling ionization 
process. 

According to the simpleman model, the electron appears in 
continuum with a vanishing momentum in the most probable 
trajectory and afterwards is driven by the laser field, absorbing 
n laser photons. The energy conservation for the ionization 
process reads: 

(Se - nieC^) + (fij - OTjC^) -flOJ- Ip, (1) 

where ee,i — c ^p^ ■ -i- .c^ and , are the energy and mass 

of the electron and the ion, respectively. The momentum con¬ 
servation provides 

Pf±+P,± = 0 (2) 

no) 

Pez+Piz = -> ( 3 ) 

c 

where , j and ,■ are the electron and the ion momentum 
components along the laser propagation direction z and the 
transverse direction, respectively. Combining Eqs. ([T]) and ([^ 
we derive 

(Be - CPez “ nieC^) + (Si - Cpiz - niiC^) = -Ip. (4) 

In the case of a common atom (m, » nie) the kinetic energy 
of the ion is small and can be neglected, e, - m,c^ =» 0. In the 
plane laser field 

K = s - cpz (5) 

is an integral of motion 0. In the simpleman model Ag - 
irigC^ for the electron because of vanishing of the electron 
kinetic momentum at the tunnel exit. Consequently, the ion 
momentum is derived 



In the simpleman model the electron transverse kinetic mo¬ 
mentum in continuum is derived from the conservation of the 
transverse canonical momentum p^j^ - eA{(p) - const. Tak¬ 
ing into account that at the ionization moment (po the electron 
kinetic momentum is vanishing, we have 


The structure of the paper is the following. In Sec. [n|the 
momentum partition is analysed using the simpleman model 
along with the energy-momentum conservation law. In Sec. 
I^the result of the strong field approximation (SFA) is pre¬ 
sented, which is followed by the discussion in Sec. IV and the 
conclusion in Sec. lYl 


[A((^o) - A((^)] , 


(7) 


where (p - coit - zjc) is the laser phase. The electron longitudi¬ 
nal momentum is derived from the conservation law 


K,i + CPe.iz = C yjpl^ 


2 2 


( 8 ) 
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which yields 


c^pl+m]/-Kl. 

2che i 


(9) 


Assuming a vanishing transverse momentum at the tunnel exit 
(Af = nieC^) yields 


Pez = — [A(<^o) - A((^)] . 

Zfilg 

(10) 

At switching off the laser field A((f>) 0: 


P.r = -eA(0o). 

(11) 

Up 


Pez = —: 

(12) 

c 


Eg - MgC^ + Up, 

(13) 


where Up - A^(^q)I(I nig) - Therefore, accord¬ 

ing to the simpleman model, Eqs. Q, (j^ and ( |12[ ), from the 
totally absorbed photon momentum during ionization 


no) 


c 



(14) 


the ion absorbs the momentum /7,j. - Ipjc and the rest of the 
photon momentum Upjc is transferred to the electron. 

The experiment of Ref im is in accordance with the above 
mentioned results of the simpleman model. However, the 
ion dynamics during the tunneling and the relativistic fea¬ 
tures of the electron under-the-barrier dynamics ifTOl fT4l [TSll 
have impact on the momentum partition in strong field ion¬ 
ization which will be discussed in the next sections. Here we 
shortly mention that when going beyond the simpleman model, 
the Lorentz force during the under-the-barrier motion induces 
a nonvanishing momentum component for the electron, pf^, 
along the laser propagation direction at the tunnel exit. Then, 
Ag ^ nigC^ - cpf^, when pf^ « m^c as usually is the case, and 
from Eq. Q one derives 

Pi,^Iplc-pf^. (15) 


In the case of exotic atoms, when the constituents masses are 
of the same order, the simpleman condition for the vanishing 
momentum at the tunnel exit concerns the relative momentum 


= (m,pf^ - mg^f'’)IM = 0, (16) 


see Table I. Then, = 0, because of Eq. (|^. 

taking into account that pf^ - (milmg)pfl, one has 

Moreover, 

A; - m,c = —(Ae - MgC ), 

(17) 

Mg 


and from Eq. (|^ one derives 



(18) 


(19) 


where M - nig + nii. The electron and ion momenta can be 
calculated from Eqs. Q and ( |1 l| l, which finally leads to the 
following expressions in the leading order of Ipl{Mc^): 


Per 

Pez 

P,r 


Piz 


-eA(0o), 




eA(4>o ), 


me 




Mc^ 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


This general expression for the final electron and ion momenta 
according to the simpleman model includes the limiting cases 
of the infinitely heavy ionic core [m, ^ oo, see Eqs. (©,([TT) 
and ( |12| l] and of the positronium atom [irig = m,]. In the 
latter case there is symmetry between the electron and positron 
dynamics stemming from a positronium atom 


Pez = Piz ~ 





(24) 


where p,- denotes the positron momentum. 

Thus, the simpleman model along with the energy- 
momenrnm conservation provides information on momentum 
partition between the electron and ion in the tunneling ioniza¬ 
tion process, given by Eqs. (201-(23 i. We underline that the 
simpleman model does not take into account the initial momen¬ 
tum of the electron at the tunnel exit, which can arise due to 
the Lorentz force Eol and due to nonadiabatic dynamics (in 
the intermediate regime between the tunneling and mutiphoton 
regimes) ll29l . Moreover, the Lorentz force effect depends 
on the ionic recoil inducing corrections of the order of OTe/m,. 
Note also that the final momentum distribution is disturbed 
also by the Coulomb focusing during the electron motion in 
the continuum M- 


In the next sections the Lorentz force effect and the impact 
of the ion recoil for the momentum partitioning between the 
ion and the electron are discussed. 


III. STRONG FIELD APPROXIMATION 


We consider strong field ionization of a simple atomic sys¬ 
tem consisting of a positively and a negatively charged particles 
(labelled as ”ion” and ’’electron”, respectively) which are ini¬ 
tially in the ground state of the bound system. The main aim 
is to study the influence of the ion motion on the electron 
tunneling dynamics and its impact on the momentum parti¬ 
tioning between the ionized electron and the atomic core in 
tunnelling ionization. The theory will be applied for the cases 
of ionization of hydrogen, muonium, muonic hydrogen, and 
for positronium. The degree of freedom both of the ion and 
the electron will be taken into account. The effect of the mag¬ 
netically induced Lorentz force on the ionization dynamics is 
included by means of nondipole treatment of the laser field in 
the weakly relativistic regime. 
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The Hamiltonian of the system under the consideration, 
consisting of an ion and an electron in a laser held reads: 



1 

> 

1 

1 

ri 1 L" 

1 

+ - 

2OT(, 

1 

< 

1 

'4) 


(25) 


where rej, qt — 1 and qe - are the radius-vectors and 
charges of the ion and the electron, respectively (henceforth 
atomic units are used). To simplify the calculations it is conve¬ 
nient to write the Hamiltonian in the relative and the center-of- 
mass (c.m.) coordinates: 



The variable transformation is shown in Table I. We use 
nondipole description in order to study the dynamic of the 
system under the inhuence of the magnetic component of the 
laser held which is responsible for the momentum transfer 
from the photons to the electron and the ion. The circularly 
polarized laser held propagating in z-direction is described by 
a vector potential: 


momentum of c.m. of the electron-ion system is Pq; |'Tp,p(f)) 
is the continuum state of the electron and ion in the laser held 
with the asymptotic c.m. momentum P and the relative mo¬ 
mentum p, neglecting Coulomb interaction. 

Similar to Ha, the transition matrix element of Eq. ( [28] l can 
be represented as 


Mn = 


X oo 

dt<T'p,p(f)| y(r)|T'o,p„(f)>, (29) 

oo 


where y(r) = qiqelr is the atomic potential. 

The continuum wave function for the electron and ion sys¬ 
tem in the laser held, after transformation to the light-time 
T - t - Zjc, where Z is the c.m. coordinate along the laser 
propagation direction, fulhls the equation 


/5.'I'p,p(t) = HT'p.pCt) , (30) 


with 


1 

2ot; 



•T7 K I 

-I- i V - <7, A T H- 

^ \ Me! 





•T7 a/ 

^ \ Mcl 


Taking into account the conservation law [P, H] - 0, the c.m. 
coordinates are factorized in the wave function: 


Ax(Za,t) = Aq COS (Jj(t-Zalc) 

Ay{Za,f} - Ao sin 0 j{t - ZaIc) (27) 


with Ao = EqIu) and a 6 {i, e]. 

The ionization transition amplitude is calculated using SEA 

II30H321: 


X oo 

df<'Tp,p(f) I Vdt) I 'To,P„(f)> 

oo 


(28) 


where 


vlw =x; 

a 


me 


p„A(z„, f) 


qlA^iZg, t) 

Inig 


describes the interaction with the laser held; |T'o,Po(f)) = 
|(j)o)e' 4 '+'i’oK jjjg initial bound state of the electron-ion sys¬ 
tem which is in the ground state |<I)o) with the energy -Ip-, the 


R = (X, L Z) = 

m, -1- m. 

mg 

r, = R-r 

M 

Zi mg 

t -= r-l- —z 

c Me 

r = (x,y,z) = -r,- 

r. = R H-r 

M 

^ -= TT^ 

c Me 

minie 

m, -1- 

V \7 V 

M ^ '■ 

Z 

T = t - 

C 

M = nij + nig 

^ e — ~r7^R + y/- 

M 

N 

II 

ml - trig 

rj - 

mjme 

p = p. + p, 

d _ d d 

dZ edr 

P = -iVjt, p = -iV, p = yuf 

m,p^ - m^p, 

“ M 

_ dl _ dr _ 


TABLE I: The variable transformation to the relative and c.m. coordi¬ 
nates and, further, to the light-time. 


T'p p(t, r, R) = exp (iP ■ R - i&t) (^(r, t). 


(31) 


with the c.m. energy £ and momentum P, after the interaction 
is switched off. 

In the weakly-relativistic regime the vector potential can be 
expanded to the hrst order in 1 /c: 


./ maZ\ ... nigZ^^ , 

A T +- A(t) +-E(r). 

\ Mcl ^ ^ Me 


(32) 


Then, neglecting the high-order terms over 1/c in Eq. (301, we 
arrive at the following equation: 




-H — [-/V + A(t)] A'(t) 1 - — 

IfleC \ Mi 

2M Me J 

When the wave function is parametrized as 

0(r, t) = exp[/p ■ r -h !T(t)z - iS (t)] , 
one obtains 

^ p \ f [P + A(t')] A'(t') , 


(33) 


(34) 


(35) 


S(t) = 


^ P . f dT'[p + zT(T' 

1 - 4 


)+A(t')]" ,(36) 


2 /^( 1 -I 7 ) 
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where ^ and 77 are defined in Table I. Finally, the continuum 
wave function for the electron-ion system in the laser field in 
the leading order of the 1 /c-expansion is 


'Pp.pCr, r, R) = Af exp (iS) , (37) 


with the normalization constant N and the action 

A 2 (t) 


S - P R-£f-i-p r-i- 


Me 


P ■ A(t) h- 


yU \ Me ! J_oo 


P ■ A(t') -h 


2 

A 2 (t') 


<38) 


We calculate the matrix elements of the amplitude between the 
ground state and the continuum as: 


Mfi = 


-iN JdtJdRjdr 


e-‘-^V(r)0o(r)e‘^'’', (39) 


assuming that the atom is at rest in the initial state Pq = 0. 

Let us first consider the momenrnm sharing between the 
electron and the ion in the simplest and transparent case when 
the atomic potential is modelled by a short-range potential. 
Later we will discuss the correction to this picture due to the 
real atomic potential. In the case of short-range potential 


where 


O-Q 




A - r ^ ^ ^ P' 

A = Ip - nto -I- ctq H- -- 1 - -— , 

Znig Znii 

, Pez + Piz , 7 /'«/ me \ 


2/7 
niee^Pei. 
pa> 


1 H- 


Me 

Pez 


Me 


Piz JLl'Hi _ !!!£ ' 


After the integration over the ion momenta, one obtains for the 
photoelectron momentum distribution 


dW 


d^Pe (277)4 J 




(47) 


Here we take into account that 


dp, 


Pez Piz 


nco — (To 


= 1 -H ^- 1 . 


Expressing the ion momentum via the electron momentum 
using the momentum conservation 5-function, we have 


V(r) = {Znj K)6{r)drr, 


(40) 


one has (p|L|<I>o) = - V^/(27r), with k = ^jZjjTp. We expand 
the last term of the action in Eq. 


expi--(l + 


(■ 


Pz + Wz 
Me 






^ exp {inojT + icr qt + intpo) , 


(41) 


where Jn(0 the Bessel functions, n is the number of absorbed 
photons, tam/jQ = PxIPy, and 


J 1 Pz+PPz'' 

j P±Ao 

(42) 

Me ) 


J 1 Pz + ^Pz^ 

42 

(43) 

Me i 

'2p 


o-Q = ( 

and obtain for the transition amplimde 

00 

Mfi = /(27r)3^V^5(P^) ^ h(Pz + 


(To na)\ 
e e 1 


X 6 


‘ r? 

, P 


2M 2 ^ 


lAOi 


inipQ 


(44) 


Then we can calculate the ionization rate 


V = 


lUee^Pei. 

pai 

1 

i + f^4 

^ mr 


1 + 11 -^ 


nii j niee Mnij 


\ no) Y me 


Pez 

\ 

( irige^ 2 p 

\ '«(■ ( 

nieC 

J 


(48) 

(49) 


A 


Zpmee^ \ myc / 2 m, j 


(50) 


where Ip = Iplniee^ and v = (nw - cro)fmee^ . 

The ionization rate given by Eq. depends on ion mass 
and takes into account the impact of the motion of the c.m. of 
the electron-ion system on the tunneling dynamics. Eor the 
hydrogen atom it is negligible as it scales with the small ratio 
melnij ss 1/1836. However, for exotic atoms this effect cannot 
be neglected. 


IV. DISCUSSION 

Let us analyse the ionization differential rate to find out the 
most probable momentum of the ionized electron and ion in 
the case of different atomic systems. We can approximately 
replace the summation over the photon number n in Eq. 
by integration and carry out the latter using the 5-function: 


dW = \Mfif- 


d^P d^p 

(Zny (27r)^ 


= \Mpf 


d^Pi d^Pe 

(2ny (27 t)^ 


(45) 


Describing the final phase space via the electron and ion mo¬ 
menta, the ionization rate reads 


dW 

d^Pid^Pe 


77 V' ^ . J noj-a-Q\ 

^ ^ n 

x5(A)72(^), (46) 


dW K 

d^Pe (27r)4 / Pe^ nie Y 

\ m;c nii j 

where we have used that 

M _ u ^ _ //L \ 

dn I ^ 7 «,c ifii j 


(51) 


( 52 ) 
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Here, the number of absorbed photons, or the parameter v, 
is determined from the energy conservation A = 0, whose 
approximate solution reads 


V 


/ + 


1 + 


Pez 


(53) 


where we have used that Ip ^ I and pi « pm^c^. Accordingly, 
the number of absorbed laser photons is 


obtain: 


Pel. 


Pez 


nigC^ 




pL 

2m„c 


(61) 

(62) 


In the latter the leading order terms with respect to y are re¬ 
tained {y - < 1 in the tunnelling regime). The ion momen¬ 

tum can be deduced from the d-functions of Eq. (|46]l. 


nieC 

OJ 


2\ 


1 -H ^- zr 


fM 

2 m,' 


1 + ii-^U^ 

mi} rUeC 


(54) 


The qualitative behaviour of the momentum distribution ac¬ 


cording to the differential ionization rate of Eq. (51 1 is deter¬ 


mined by the Bessel function. In the tunneling regime, when 
t/p/« » 1 and /p » 1, one has » 1, therefore, for the fur¬ 
ther analysis we will use the asymptotics of the Bessel function 
n ~ ^ ^ oo ||34l : 


Jl(0 


2n^2(n-0^ 


exp 


4 V2 (n - 1 

^ j ■ 


(55) 


Pi^ - VnigC^ - Pez ^ 
Pi± = -Pe±- 


and using Eqs. (491,(61 1 , and (621: 


Piz = 


2^ 

3c 


(63) 

(64) 


(65) 


Comparing Eq. (65 1 with our qualitative discussion in Sec. II, 
see Eq. ( [T5] l, we can conclude that the electron momentum at 

' ' /O') 

the tunnel exit is - Ip 13c, which is the reason of variation 
the ion momenrnm from the Ipjc value. This result coincides 
with the predictions of Refs. ifTOlfT?! . 


The peak of the momenrnm distribution corresponds to the 
minimum of the expression (F = (n - Cp' achievable at 
n ~ where n and ^ are given by Eqs. ([54)i and (48 i. 


B. Exotic atoms 


A. Hydrogen atom 


Eirst, let us consider the simplest limit of infinitely heavy 
ion m,' —> oo. In this case. 


V = 

n = 


x^ 

— + 
2 

nieC^ 

a> 

nieC^ 


to 



^R±(l +Pz), 

. ec^+Pz) 

y + —;- 


(56) 

(57) 

(58) 


with;><:^/2 s ip,p±_ = pex/mgC and p^ = pej/mgC, which yields 


T(pi_,Pz) 


nieC^ +pI+pI + (^^ - 2p±^)(l + Pz)]^ 
w j 8fpi(lH-pz) 


(59) 

The conditions for F = min, &F/dp_i_ - 0 and 5F jdp^ - 0, 
read, respectively: 


In the case of exotic atoms, such as positronium (nielmi - 
I), muonic hydrogen atom {niglmi 0.II26), and muoninm 
{nielmi 1/207), the masses of constituents are comparable 
and, therefore we have to use the general expressions for the 
parameters n and given by Eqs. (|48|,(|4^, and (54 1 . 

Eirst, we hnd an approximate solution of Eq. (49 1 in pertur¬ 
bation with respect to the parameter 


V 


v( 2 ) + v(4) , 


( 66 ) 


where ~ assuming that p± ~ ^ and p^ ~ 

v( 2 ) ^ d(i+^)+^, (67) 

y(.4) ^ (1 + ^ ( 68 ) 


where p s nielnii. The parameters of the Bessel functions up 
to the order of are: 


6 [pi. - ^(1 + Pz)] P± ^ + pI + pI + if - 2p±^)(l -H pj, 

3(l+p,)[2pz + ^2_2^^^j^ (60) 

^x^+pI_+pI + (f - 2p^^)(l -H pj. 

Solving the latter in perturbation with respect to ^ and taking 
into account that pi_ ~ ^ and p^ ~ as well as x ~ we 


to 


nigC 


-(1 H-p)^Pi 


1 -H(l -p)pz + 


\ + p 


./ 2 ) 


, (69) 


1 


ff 



( 70 ) 
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The condition d'F I dp i_ - 0 in this case yields: 


3/?±(l +p) 


1 +(\-p)p, + 


1 + ^ 


{(1 +fp^)pi 


= y(2) 




1 + 


1 + (1 - p)p^ + + piiP 

I + p 




- (1 +P)^Pi 


1 +/i + (l + -^y® 

I + p 


1 + (1 - p)p^ + + piiP 

i + p 


which in the leading order reads 


(71) 


at the tunnel exit. The longitudinal component of the ion 
momentum, then, is 


Piz 


Pi. me 
InipC nii 


3c 


(l - — 

\ 2M) 


ill) 


The electron and ion longitudinal momentum Eqs. (761 and 


( [77| are different from the prediction of the simpleman model 
Eqs. ( |21| l and (23 i. It is due to the nonvanishing electron-ion 
relative momentum at the tunneling exit which depends on the 
ionization energy Ip as well as on the mass ratio iriglM. The 
latter factor describes the role of the ion motion during the 
tunneling process. 

Now, from Eqs. ( [73| ), (761, and \11\ , we are able to eval¬ 
uate the most probable momentum for ionization of exotic 
atoms, taking into account the effect of the ion motion on the 
ionization dynamics. 


3p±(p± - f) =-T-+ 


2(1+p) 


(72) 


1. Muonium 


The solution of the latter provides us the most probable trans¬ 
verse momentum: 


Pel. = nieC^ 


6 1 + 


(73) 


The derived transverse momentum component contains a nona- 
diabatic correction, the term ~ in Eq. (73 i, which is absent 


in our simpleman estimation via Eqs. (201 and which is dis- 


mrbed by the ion recoil (see the term ;««/;«,). 

The second condition of the maximal probability = 

Ois 


In the case of muonium (electron and antimuon) the most 
probable momentum of the electron is 


Pe± - (l + 0 . 1667 ^) , 

Pel. , ci 

Pez ~ 0.335 . 

InieC C 

The momentum of the antimuon is 

Ppi. ~ OTeCf ( 1 -H 0.166y^) , 
Ip 


Ppz 


+ 0 . 665 '" 


414ot»c 


(78) 

(79) 

(80) 
(81) 


l+(l-p)p,+ ^P‘ 

1 + p 


(1 - fP)^ + Pzd +P)- 


-(i-~p^)^Pi]^(i-~p)U^Hi + 


-H(l -H/i)[l +(\-p)p^] j 


(1 +P)^P± 


1 +)! + (1 -p)Pz + 




,/ 2 ) 


1 + p 


2. Muonic hydrogen 

In the case of a muonic hydrogen atom (muon and proton) 
the most probable momenrnm of the muon is 


(74) 


In the leading order, the Eq. ( [75] > is simplified: 

,y2 




-(1 -p)^Pi 


Ppi. ^ m<,c^ (1 H-0.14987) , 

(82) 

Ppz ~ T +0.367 , 

2mc 

(83) 

While for the proton they are 


Ppi. ~ (1 H- 0 . 14987 ^) , 

(84) 

Ppz - -|^+0.633^, 
lom^c c 

(85) 


= (1-/1) 


/2i 


^2 

- + - +- ^Pl. 

2 2(1+p) 2 


(75) 


the solution of which provides us the most probable longitudi¬ 
nal momentum: 


pIi. 

Pez ~ ^ 

InieC 


3c\ Ml 


3. Positronium 

In the case of a positronium atom the most probable momen¬ 
tum of the photoelectron is 


(76) 


The ion momentum is derived from Eqs. (631 and (641. The 
second term corresponds to the momentum of the electron 


Pel. 

= m.c^(l + ^), 

(86) 

Pez 

_ ^ Ip 

2meC 2c 

(87) 
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The positron momentum components are the same by the ab¬ 
solute value (the transverse momentum is opposite). 


C. The role of the Coulomhic atomic potential 


In the discussion above, we assumed a short-range atomic 


potential, Eq. (401. Now we examine how the SEA calculations 
are modified when the exact Coulombic atomic potential is em¬ 
ployed. In this case the matrix element (p|y|Oo) = - ^l{2n) 
should be replaced by 


<p|V|<to) 


4 


K^+p\ + [p^ -H/?(t)] 


4 

+ p^ 


1 - 


2pJ3(t) 


+ p^ 


( 88 ) 


where a - jti and 


j0(t) = Pq+Px sin(wT + (fio) 


Pa 

Pi 


^ niecf 


nil 


m„ , 

1-— KEx 

mi 


In Eq. ( 881 we have expanded the expression with respect 


to j8(t) ~ 1/c. Then rather than Eq. ( |46| l, we will have the 
following differential ionization rate 


dW 




2 3 

a 


n 4—i {k^ + p^p 


b(PJ(5 P,+ 




X (5(A) JliO 


1 - 


^Pz I ^ 


+ p^ 


c 


(89) 


The ionization differential rate integrated over the ion momenta 
reads 


dW 

d^Pe 


X 


4 

^1 


JliO 


Pez 


Ml 




( k ^ + p2)2 


Adj I n 

-T^ + -A 

+ P^\ Z, 


(90) 


where p^ - pl^ + i^p^z - rn^cvIM^ and p^ - 
The parameters v, n are determined by Eqs. 


•HI. 


Pez - m^cvIM. 
( |4^ , ( [39| and 


The qualitative behaviour of the ionization rate is illustrated 
in Fig. [T] on the example of the positronium ionization. The 
accounting of the exact atomic potential corrects only the value 
for the ionization rate, but the position of the peak of the 
momentum distribution is determined by the Bessel function. 
Therefore, one can rely on the conclusions on the momentum 
sharing between the ion and the electron presented in this sec¬ 
tion above. Figure 1 illustrates that our analytical expression 


o Eq. ( |87|l [w hich is the particular nig - nii case of the gen¬ 
eral Eq. (|76|l] provides the correct value for the longitudinal 
component of the electron momentum at the maximum of the 
distribution. 

KCdW/rf^p) 



FIG. 1: (color online) The electron distribution function over the 
momentum along the laser propagation direction = PezlinieC^^ 12) 
in the case of a positronium ionization: (solid) via Eq. ( |90| > using 
Coulomb atomic potential; (dashed) via Eq. 0 (multipled by a 
factor of 100) when atomic potential is modelled by a short range 
potential. The transverse momentum is evaluated at the maximum 
of the momentum distribution = ingC^il + 7^/12) according to 
Eq. l |86^ , Pz,max = (Pexlnic^)^ + 7^/2 corresponds to the longitudinal 
component of the electron momentum at the maximum of the distri¬ 
bution according to Eq. ( |87[ l. The laser parameters are EojEa = 0.2, 
oj = 0.05 (7 = 0.5, ^ = 0.0172). 


V. CONCLUSION 

We have investigated the momentum partition between the 
constituents of exotic atoms during strong field tunneling ion¬ 
ization. The momentum distribution is deviated from the pre¬ 
diction of the simpleman model. One reason for the deviation 
is that the electron appears in the continuum with nonvanishing 
momentum along the laser propagation direction which is due 
to the effect of the magnetically induced Lorentz force during 
the under-the-barrier dynamics and due to nonadiabatic effects. 
The second reason is the impact of the recoil of the atomic core 
on the tunneling dynamics and, therefore, on the momentum 
shift of the electron (muon) along the laser propagation direc¬ 
tion. The second factor is negligible for common atoms but 
significant in the case of exotic atoms such as muonic hydrogen 
and positronium. 
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